Renormalization of chiral two-pion exchange NN interactions. 
Momentum vs. coordinate space. 
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The renormalization of the chiral np interaction in the ^So channel to N3L0 in Weinberg counting 
for the long distance potential with one single momentum and energy independent counterterm 
is carried out. This renormalization scheme yields finite and unique results and is free of short 
distance off-shell ambiguities. We observe good convergence in the entire elastic range below pion 
production threshold and find that there are some small physical effects missing in the purely pionic 
chiral NN potential with or without inclusion of explicit A degrees of freedom. We also study 
the renormalizability of the standard Weinberg counting at NLO and N2L0 when a momentum 
dependent polynomial counterterm is included. Our numerical results suggest that the inclusion of 
this counterterm does not yield a convergent amplitude (at NLO and N2LO). 
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I. INTRODUCTION 



The modern Effective Field Theory (EFT) analysis of 
the NN interaction using chiral symmetry as a constraint 
has a recent but prolific history |3] (for comprehen- 
sive reviews see e.g. Ref. [3, 0, Q)- Most theoretical 
setups are invariably based on a perturbative determina- 
tion of the chiral potential [1, 0,11, i, [13, [HI, [ll [ll [11 
and the subsequent solution of the scattering prob- 
lem [H [H H [H H mi, m m, H p. Actu- 
ally, the theory encounters many problems in the low 
partial waves and in particular in the s— waves (see how- 
ever Ref. [IHl for a more optimistic view). Indeed, there 
is at present an ongoing debate on how an EFT pro- 
gram should be sensibly implemented within the NN con- 
text and so far no consensus has been achieved (see e.g. 
[23, m, H [S^])- The discussion is concerned with the 
issue of renormalization vs. finite cutoffs, a priori (power 
counting) vs. a posteriori error estimates or the appli- 
cability of perturbation theory both on a purely short 
distance theory or around some non-perturbative dis- 
torted waves. At the moment, it seems fairly clear that 
an EFT scheme with a cutoff-independent and system- 
atic perturbative power counting for the S-matrix (the 
so-called KSW counting) fails [sil, [sl, [H, [H. On the 
other hand, the original EFT inspired scheme P, (the 
so-called Weinberg counting) has recently been shown 
to produce many results which turn out to be strongly 



cutoff dependent [53, [1^ and hence to be incompatible 
with renormalizability ^. Thus, some acceptable com- 
promise must be made. Actually, the chosen approach 
between this dichotomy depends strongly on the pursued 
goals and it is fair to say that any choice has both ad- 
vantages and disadvantages. In any case, the reason for 
both failures can be traced back to the nature of the long- 
distance chiral potentials; while pion-exchange potentials 
fall off exponentially at long distances they include strong 
power-law singularities at short distances. Those sin- 
gularities become significant already at distances com- 
parable with the smallest de Broglie wavelength probed 
in NN scattering below pion production threshold. Ob- 
viously, any development of NN interactions based on 
chiral dynamics will presumably require a deeper under- 
standing and proper interpretation of the peculiarities of 
these highly singular chiral potentials. Although singu- 
lar potentials where first analyzed many years ago (36[ 
(for an early review see e.g. [33| and for a more updated 
view within an EFT context see [sl]), their short dis- 
tance singular character within the NN interaction has 
seriously been faced more recently within a renormal- 
ization context for the one-pion exchange (OPE) poten- 
tial 2a, M El, m, Ii3|j4l&5^ and the two-pion 
exchange (TPE) potentialillllJljT 



In the np scattering problem the ^6*0 channel is very 
special since the scattering length is unnaturally large 
as compared to the range of the strong interaction, 
ao = -23.74(2) fm > I/m^ = 1.4 fm. In fact, even 
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^ In this paper we refer to renormalizability in the sense of param- 
eterizing the short distance physics by a potential which matrix 
elements in momentum space are a polynomial in the momenta. 
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at zero energy the wave function probes relatively short 
distance components of the chiral potential [35| ^. Higher 
energies become even more sensitive to short distance in- 
teractions. Consequently, this channel looks like an ideal 
place to learn about the size of the most relevant short 
range corrections to the NN force in the elastic scatter- 
ing region. Actually, in the ^5*0 channel, most EFT in- 
spired schemes yield at leading order (LO) (which con- 
sists of OPE plus a nonderivative counterterm) an al- 
most constant phase of about 75° around k = 250 MeV, 
and an effective range of = 1.44fm. On the other 
hand, most determinations from Partial Wave Analy- 
sis [43, [l^l and high quality potential models [5l| yield an 
almost vanishing phase at this center-of-mass (CM) mo- 
mentum while the experimental effective range is about 



twice the OPE value, 



2.77(5)fm. 3 It is quite 



unbelievable that such large changes can be reliably ac- 
commodated by perturbation theory starting from this 
LO result despite previous unsuccessful attempts treat- 
ing OPE and TPE perturbatively [Sll [H 113, [H . Ac- 
tually, for the singlet channel case, short distance com- 
ponents are enhanced due to the large value of the scat- 
tering length and the weakness of the OPE interaction 
in this channel. This is why TPE contributions have 
been treated with more success in a non-perturbative 

fashion I, M, M [13, M M M iH, m, H, Hi HI- 

Despite phenomenological agreement with the data, the 
inclusion of finite cutoffs suggests that there might be 
some regulator dependence in those calculations. 

In a series of recent papers [H, 13 El], two of us 
(M.P.V. and E.R.A.) have proposed not only to iterate 
but also to rcnormalize to all orders the NN chiral po- 
tential within a long distance expansion. By allowing the 
minimal number of counterterms to yield a finite result, 
long-distance regulator-independent correlations are es- 
tablished. In practice, the potential must be computed 
within some power counting scheme. While the potential 
is used within Weinberg's power counting to LO, next- 
to-leading order (NLO), and next-to-next-to- leading or- 
der (N2L0), we only allow for those counterterms which 
yield a finite and unique scattering amplitude. In the ^Sq 
np channel a single energy and momentum independent 
counterterm Co is considered which is determined by ad- 



This can be best seen by means of the effective range formula 

2" 



Jo 



dr 



ao 



where no(^) is the zero energy wave function, fulfilling the asymp- 
totic condition uo{r) — > 1— r/oo- Most of the integrand is located 
in the region around r = Ifm which is in between OPE and TPE 
ranges. Moreover, the low energy theorem of Ref. [ssl allows to 
write ro = A-i- B/ao + C/ag which in the extreme limit ao — ► 00 
yields rg — > A. Numerically it is found that A is far more de- 
pendent than B and C when being evaluated at LO, NLO and 
NNLO. 

In fact, these high quality potential models yield slightly smaller 
values, ro ~ 2.67fm. 



justing the physical scattering length. In the derivation 
of this result the mathematical requirements of complete- 
ness and self-adjointness for the renormalized quantum 
mechanical problem for a local chiral potential play a de- 
cisive role. This surprising result is in contrast to the 
standard Weinberg counting where an additional coun- 
terterm C2 is included already at NLO. This C2 coun- 
terterm could, in fact, be determined by fitting the ex- 
perimental value of the effective range; the physics of 
C2 is to provide a short distance contribution to the ef- 
fective range in addition to the contribution from the 
known long distance chiral potential. Within this con- 
text it is remarkable that according to Ref. [3^, where 
such a short distance contribution vanishes (or equiva- 
lently C2 — when the cut-off is removed), rather ac- 
curate values are predicted yielding r^'"*^ — 2.29 fm and 
^N2LO _ 2.86 fm after renormalization. This latter value 
is less than 3% larger than the experimentally accepted 
value and it suggests that most of the effective range 
is saturated by N2L0 TPE contributions and calls for 
pinning down the remaining discrepancy. This trend to 
convergence and agreement is also shared by higher order 
slope parameters in the effective range expansion with- 
out strong need of specific counterterms although there is 
still room for improvement. Motivated by this encourag- 
ing result, one goal of the present paper is to analyze the 
size of the next-to- next-to-next-to- leading order (N3L0) 
corrections to the results found in [ssj . 

The calculations in Ref. [H, Hi, [H, [s^ exploit explic- 
itly the local character of the chiral potential by conduct- 
ing the calculations in coordinate space which makes the 
analysis more transparent. Many results, in particular 
the conditions under which a renormalized limit exists, 
can be established a priori analytically. Moreover, the 
highly oscillatory character of wave functions at short 
distances is treated numerically using efficient adaptive- 
step differential equations techniques. This situation con- 
trasts with momentum space calculations where, with the 
exception of the pion-less theory, there is a paucity of 
analytical results, and one must mostly rely on numeri- 
cal methods. Moreover, the existence of a renormalized 
limit is not obvious a priori and one may have to resort 
to some trial and error to search for counterterms. Fi- 
nally, renormalization conditions are most naturally for- 
mulated at zero energy for which the momentum space 
treatment may be challenging, at times. Of course, be- 
sides these technical issues, there is no fundamental dif- 
ference between proceeding in momentum or coordinate 
space, particularly after renormalization, provided the 
same renormalization conditions are specified, since dis- 
parate regulators stemming from either space are effec- 
tively removed. Indeed, we will check agreement for the 
phase-shifts determined in different spaces whenever such 
a comparison becomes possible. This equivalence is in it- 
self a good motivation for renormalization. 

However, at N3L0 some unavoidable non-localities ap- 
pear in the chiral long distance potential. Although they 
could be treated in configuration space, we adopt here a 



3 



momentum space treatment. This will also allow us to 
answer an intriguing question which was left open in the 
coordinate space analysis of previous works [s^, [13, \^ , 
namely, the role played by the conventional momentum- 
polynomial representations of the short distance interac- 
tion used in most calculations 

IHj H^j HE] in the renormalization problem. More specifi- 
cally, Ref. [IHl showed that taking C2 = was consistent, 
and a regularization scheme exists where a fixed C2 was 
irrelevant, but could not discriminate whether C2 ^ 
was inconsistent as far as it was readjusted to the ef- 
fective range parameter for any cut-off value. ^ In this 
regard, the present paper yields a definite answer mak- 
ing the surp rising agreement of the effective range found 
in Ref. [35] an inevitable mathematical consequence of 
renormalization. 

The paper is organized as follows. The problem is 
stated in Sec. |TT] where a general overview of the renor- 
malization problem is given both in momentum as well 
as coordinate space. In Sec. Ill B[ we particularize the 
momentum space formulation of the scattering problem 
with counterterms for the ^Sq channel within a sharp 
three-momentum cutoff scheme. Likewise, in Sec. Ill CI 
we proceed similarly in the coordinate space formula- 
tion within a boundary condition regularization with a 
short distance cutoff Vc- In Sec. IIIII we discuss some 
features of both coordinate and momentum space for- 
mulations in the pion-less theory and try to connect the 
high momentum cutoff A with the short distance radial 
cutoff re- This allows a one-to-one mapping of coun- 
terterms in both spaces which will prove useful later on 
in the pion-fuU theory. The identification between the 
sharp momentum cutoff and the short distance radius 
found in the pion-less theory is discussed further in Ap- 
pendix 1^ in the presence of a long distance potential in 
the light of the Nyquist theorem. In Sec. IIVI we come to 
the central discussion on N3L0 corrections to the phase 
shifts when the scattering amplitude is renormalized with 
only one short distance counterterm. A wider perspec- 
tive is achieved by further considering the role of explicit 
A-excitations in intermediate states and the subsequent 
one-counterterm renormalization of the scattering ampli- 
tude. We also discuss the role of three pion exchange as 
well as how the results depend on the renormalization 
scheme used to compute the potential based on cut-off 
independent counterterms. The mathematical justifica- 
tion for using just one counterterm is provided in Sec. |V] 
where the standard Weinberg scheme is pursued both in 



* We mean of course the case when the cut-off is being removed. 
The essential issue is whether or not one can fix by a short dis- 
tance potential which is a polynomial in the momenta the effec- 
tive range independently on the potential and remove the cut-off 
at the same time. Of course, the very definition of the poten- 
tial is ambiguous and requires a specific choice on the polynomial 
parts. Technically, we find that any fixed, cut-off indendependent 
C2, becomes irrelevant in the limit A ^ 00 (see below). 



momentum and coordinate space at NLO and N2L0 and 
shown to potentially have some problems. Finally, in 
Sec. IVH we summarize our main points. 

II. THE RENORMALIZATION PROBLEM 

A. General overview and main results 

Let us define the scope and goals of the present work. 
The standard non-perturbative formulation of the renor- 
malization problem starts with an effective Lagrangian 
or Hamiltonian (see e.g. [E and [3, Q and references 
therein) , from which a certain set of irreducible Feynman 
diagrams (usually up to a certain order) is calculated. 
These irreducible diagrams are defined to represent a po- 
tential V. The potential is then inserted into a scattering 
equation where it is iterated infinitely many times or, in 
other words, re-summed non-perturbatively. In the CM 
frame, where the np kinetic energy is given hy E = p^/M, 
with M = 2^„p = 2M„Mp/{Mp + M„), the scattering 
process is governed by the Lippmann-Schwinger equation 

T^V + VGoT, (1) 

with V the potential operator and Go = {E — Hq)~'^ the 
resolvent of the free Hamiltonian. The outgoing bound- 
ary condition corresponds to E ^ E + iO+. Using the 
normalization {x\k) — e''^'^/(27r)'^/^ one has 

(2) 

Here A means a generic regulator and represents the scale 
below which all physical effects are taken into account 
explicitly. The degrees of freedom which are above A 
are taken into account implicitly by including a suitable 
cutoff dependence in the potential. The precise equation 
governing this cutoff dependence was described and stud- 
ied in some detail in Ref. [s^ with particular emphasis 
on infrared fixed points. We will analyze the cutoff de- 
pendence below focusing on the ultraviolet aspects of the 
interaction. 

Motivated by the low energy nature of the effective 
theory, the potential is usually separated into short and 
long distance components in an additive form 

{k'\V\k)^Vs{k',k) + VL{k',k), (3) 

where the long distance contribution is usually given by 
successive pion exchanges 

VL{k', k) = Vi^k', k) + l/2x(fc', fc) + . . . , (4) 

and the short distance component is characterized by a 
power series expansion in momentum 

Vs{k', k)^Go + Gik- k' + C2(P + /?') -I- . . . , (5) 
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where for simplicity we assume a spin singlet channel ^. 
Note that Cq and C2 contribute to s-waves while Ci con- 
tributes to p-waves, and so on. One should face the fact 
that, although the decomposition given by Eq. ^ is per- 
turbatively motivated and seems quite natural, the ad- 
ditivity between short and long range forces is actually 
an assumption which has important consequences, as we 
will see. ^ 

In principle, for a given regularization scheme charac- 
terized by a cutoff A, the counterterms Co, Ci, C2 and so 
on are determined by fixing some observables. One nat- 
urally expects that the number of renormalization con- 
ditions coincides with the number of counterterms in a 
way that all renormalization conditions are fully uncor- 
related. The statement of UV-renormalizability is that 
such a procedure becomes always possible when the cut- 
off A is removed by taking the limit A — > 00. This may 
not be the case as there may appear redundant contribu- 
tions (sec the discussion below in Sec. IIIip meaning that 
one counterterm or counterterm combination can take 
any value. Another possible situation is just the opposite; 
one may want to impose more renormalization conditions 
than possible. In this case some counterterms or coun- 
terterm combinations are forbidden. Rather than being 
intricate mathematical pastimes, these features have al- 
ready been investigated recently [13, [1^ for large cutoffs 
casting some doubt on the regulator independence of the 
original proposal [1,, ^ . 

Even if one admits generically Eq. ^ as well as Eq. (jH) 
and Eq. (O, it is not obvious how many terms should be 
considered and whether there is a clear way of defining a 
convergence criterium or identifying a convergence pat- 
tern. It is fairly clear that, on physical grounds, one 
should consider an expansion of the potential that starts 
at long distance and decreases in range as the number 
of exchanged particles increases. Note, however, that, 
while OPE is well-defined, the general form of the TPE 
potential is not uniquely determined. Within this con- 
text, one of the main attractive features of the EFT 
approach has been the definition of a power counting 
scheme which provides a hierarchy and an a priori cor- 
relation between long and short range physics. In the 
present paper, we will assunrc Weinberg power count- 
ing [i, m, U, i, M, M, III, liJ, El (see below for a more 
precise definition), where the long distance potential is 
determined in a dimensional power expansion. 

The long distance component Vl is obtained by par- 
ticle exchanges and, in some simple cases, depends only 



^ More detailed expressions including spin triplet channels can be 
looked up e.g. in Ref. 4j. 

® Specifically, Eq. JSjl does not foresee for instance terms of the 
form y(p',p)C(A), i.e. terms which are not polynomial but in- 
fluence the renormalization process. Of course, once additivity 
is relaxed there are many possible representations, in particular 
for the short distance components. 



on the momentum transfer ^. In such a case, if we for- 
mally take a Fourier transformation of the long distance 
potential 



Vl{x) = 



d^p 
(2^ 



VLip)e''' 



(6) 



and take the limit A ^ 00 one has the standard 
Schrodinger equation in coordinate space, 



1 

M 



(7) 



where the coordinate space potential is 

V{x) = Vl{x) + Co6^^\x)+CiV6^^\x)V 



C2 



V^5'^^^x)+6^'\x)V' 



(8) 



The whole discussion, which has been carried on for 
years now, concerns the precise meaning of these delta 
and derivatives of delta interactions, particularly when 
a long distance potential is added to the short distance 
one. A crucial finding of the present paper based on a di- 
rect analysis in momentum space is that non-perturbative 
renormalization imposes restrictions on the number of 
terms and form of the short distance potential which 
depend also on the particular long distance potential. 
Some of these restrictions were discussed in previous 
works |27i . [ssl |44| . Remarkably these new renormaliz- 
ability restrictions were conjectured in coordinate space 
in Ref. [H, IH, HI for the ^5*0 channel based on self- 
adjointness and completeness of states and apply to the 
TPE potential; a single Cq 7^ counterterm is allowed 
while two counterterms, Co 7^ and C2 7^ 0, are forbid- 
den 



B. Momentum space formulation 

In the ^So channel the scattering process is governed 
by the Lippmann-Schwinger equation 



r(fc', k) = V{k\ k) + / dqV{k', q) 



p^ — q^ + iO 



-r(q,fc) 



(9) 



where r(fc', k) and V(k' , k) are the scattering amplitude 
and the potential matrix elements, respectively, between 
off-shell momentum states k and k' in that channel and 
the sharp three-momentum cutoff A represents the scale 
below which all physical effects are taken into account 



^ This assumption will be relaxed immediately below when dis- 
cussing the Weinberg counting in the ^Sq channel. 

* Again, we mean cut-off dependent counterterms designed to fit 
physical observables. 
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explicitly. From the on-shell scattering amplitude the 
phase shift can be readily obtained 



T{p,p) = - 



nMp 



e'^ sin 6{p) 



(10) 



The short range character of the nuclear force implies 
that at low energies one has the effective range expansion 
(ERE) 



pcotS{p) = -—-Re 
Mtt 



1 



1 

ao 



1 



T{p,p) 

,2 



+ rop^+„2P +•■• (11) 



where ao is the scattering length, vq the effective range 
and V2, vs etc. are slope parameters. 

In the ^Sq channel, the potential is decomposed as the 
sum of short and long range pieces 



Vik',k)^Vsik',k) + VLik',k). 



(12) 



In the standard Weinberg counting, the short distance 
contribution is written as follows 



Vsik',k) = Co(A) + (fc2 + fc'")C2(A) 
+ Cl(A)fc2fc'' + C4(A)(fc4 + 



.(13) 



where the counting is related to the order of the momen- 
tum which appears explicitly. The long distance compo- 
nent of the potential is taken to be the sum of explicit 
pion exchanges 



(14) 



where [l| 



Vl. = 


* Itt 






^ylt' + 


V2. = 








\- ... 


V3^ = 


* 37r 


^ ... 







(15) 

using dimensional power counting. Ideally, one should 
determine the physically relevant long range regulator in- 
dependent correlations, i.e., long distance effects of simi- 
lar range. This would amount to consider all mr exchange 
effects on the same footing, since they yield a long dis- 
tance suppression ~ g-nm^r jjjodulo power corrections. 
At present, the only way how these long distance poten- 
tials can be systematically computed is by dimensional 
power counting in perturbation theory, as represented 
schematically in Eqs. and p^ . 

In the standard Weinberg counting one has 



Vlo 
Vnlo 

Vn2LO 
Vn3LO 



(0) 



* Itt 

(2) 



Vlo + K 

Vnlo + yii^ ^ 

(4) 



^ Itt 



(2) 



2. 



(3) 



Vn 



2LO 



* Itt 



V, 



(4) 



271- 



(4) 



(16) 



Note that this counting involves both unknown short- 
distance physics and chiral long-distance physics in a un- 
correlated way. Note also that there is no first order con- 
tribution and that there is no third order contribution to 
the short distance potential. 

Regarding Eq. one should stress that the separa- 
tion between long and short range contributions to the 
potential is not unique. In fact, there is a polynomial am- 
biguity in the long range part which can freely be trans- 
ferred to the short distance contribution. However, the 
non-polynomial part is unambiguous as it is directly re- 
lated to the left hand cut of the partial wave amplitude 
which for mr exchange is located at p = inmT^/2 but 
presumably becomes incomplete for \p\ > mp/2. ^ 



C. Coordinate space formulation 

In coordinate space, the problem in the ^Sq channel 
is formulated as follows [13, H^l • Assuming a local long 
distance potential Vl (r) one has to solve the Schrodinger 
equation 



Up{r) + UL{r)up{r) ^p^Up{r) , 



r > Tr 



(17) 



where Ul{t) — ^finpVLir) is the reduced potential (in 
fact, the Fourier transformation of VL{q) ) and Up{r) the 
reduced wave function for an s-wave state. Here Vc is 
the short distance cutoff and the reduced wave function 
is subject to the boundary condition at r = Tc and the 
standard long distance free particle behaviour 



Up{rc) 
Up{r) 



pcot Ss{p) , 

sin(pr + S{p)) 
sin 6 (p) 



(18) 
(19) 



where Ss (p) is the short distance phase-shift encoding the 
physics for r < Vc. In the case of a vanishing long range 
potential UL{r) — the phase shift is given by 6s{p,rc). 
On the other hand, if we take 5s{p) = we get a standard 
problem with a hard core boundary condition, Up{rc) = 
which for Tc — )■ becomes the standard regular solution 
at the origin. At low energies both the full phase-shift 
5{p) and the short distance phase-shift 5s{p) can be de- 
scribed by some low energy approximation, like e.g., an 
effective range expansion, 



pcot 5s{p) 
p cot 5{p) 



1 



1 



"0,5 2 
1 1 2 

+ T^rop 



(20) 
(21) 



^ The best way to recognize the ambiguity is in terms of the spec- 
tral function representation of the potential [7| , where the sub- 
traction constants can be fixed arbitrarily. In coordinate space 
the non-ambiguous part corresponds to the potential y(r) for 
any non- vanishing radius (see e.g. the discussion in Ref. |35| . |48|| .) 
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where ao,s is the short range scattering length, tq.s the 
short range effective range, and ao and rg the full ones 
If we also make an expansion at low energies of the re- 
duced wave function 



Up{r) = lio(r) +p^M2(r) + 
we get the hierarchy of equations 



-u'^{r) + U{r)u„{r) 
o^o,su'o{rc) + Mo (re) 
uo{r) 



0, 
0, 

r 

1 



(22) 



(23) 



and 



- (r) + U{r)u2{r) 
ao,su'2{rc) + U2{rc) 



■^rQ,sao.suo{rc) . 



(24) 



(r^ - 3aor + 3aoro) , 



and so on. The standard way to proceed would be to 
integrate the equations for uo{r), U2{r), etc. from infinity 
downwards, with a known value of ao, using Eq. (|23p to 
obtain ao,s and then one can use Eq. fTT)) together with 
Eq. (fT9| and Eq. pO|) to compute S{k) for any energy with 
a given truncated boundary condition. This procedure 
provides by construction the low energy parameters we 
started with and takes into account that the long range 
potential determines the form of the wave function at 
long distances. The only parameter in the procedure is 
the short distance radius Tc, which is eventually removed 
by taking the limit Tc ^ 0. 

One should mention at this point that the coordinate 
space is particularly suited for the case of local long dis- 
tance potentials, but the renormalization with an arbi- 
trary number of countertersm requires an energy depen- 
dent but real boundary condition at short distances which 
eventually violates self-adjointness. On the other hand, 
the momentum space formulation allows the discussion 
of nonlocal long distance potentials and the renormaliza- 
tion is done in terms of a momentum dependent short 
distance polynomial potential. Although this looks like a 
self-adjoint problem, we will see that in this formulation 
the counterterms may in fact become complex. 



III. THE RENORMALIZATION PROBLEM FOR 
THE PION-LESS THEORY 

The renormalization of the pion-less theory, i.e., a set 
of pure contact interactions, has been treated wi th g reat 
detail in the literature [ll [11, [11, [13, [11 [11 ifl [Ml, 



|62|, [63|, [6J, [65|, [66| although without much considera- 
tion on how this problem might be embedded into the 
wider and certainly more realistic situation where the 
finite range and short distance singular chiral NN po- 
tentials are present. In fact, much of the understanding 
of non-perturbative renormalization within the modern 
NN context has been tailored after those and further 
studies based on the non-singular OPE singlet ^5*0 po- 
tential [H, [13, [1^, [1^ plus the standard perturbative 
experience. In previous [H, [H, [H, [13, [1^ and in the 
present work, we pursue exactly the opposite goal: we 
will only consider renormalization procedures which can 
directly be implemented in the presence of long distance 
potentials since, after all, contact NN interactions are 
always assumed to approximate truly finite range inter- 
actions in the long wavelength limit. Thus, it is useful to 
review here those developments with an eye on the new 
ingredients which appear in the non-perturbative renor- 
malization of singular pion exchange potentials as ana- 
lyzed in later sections. In addition, the deduced running 
of the counterterms in the contact theory in the infrared 
domain serves as a useful starting point when the long 
distance pion exchange potential is switched on. Finally, 
we will also discuss the size of finite cutoff corrections to 
the renormalized result depending both on the particular 
regularization as well as the corresponding representation 
of the short distance physics. 



A. Momentum space 

Although the previously described momentum space 
framework has extensively been used in the pas t to de- 
scribe successfully the data 0, [H, [H, [H, [HI, [H, [H [13, 
[25j | with a finite cutoff A it is worth emphasizing some 
puzzling features regarding the off-shell ambiguities of 
the short distance potential when finite range corrections, 
encoded in the C2, C4 etc. counterterms, are included. 

In momentum space, the pion-less theory corresponds 
to taking VL{k',k) = 0. In such a case the Lippmann- 
Schwinger equation reduces to a simple algebraic equa- 
tion [58l[63 |. At very small values of the cutoff A < m^/2, 
the long range part of the potential may be neglected 
since they scale with powers of momentum and a sim- 
ple contact theory of the form of Eq. (fT^]) may be used. 
For instance, when Vs{k',k) = Co(A), the Lippmann- 
Schwinger equation (LSE) may be directly solved. Using 
the basic integral 



P- 



dqq^ . .Trp P A + p 

- (7^ + zO+ 2 2 A — p 



This is not the only possible short distance representation [52 
See Sec. IIII Bl for a further discussion on this. 



The result for a different momentum cutoff scheme such as 
V{k' , k) — > g{k' , A)V{k, k')g{k, A) corresponds to making the re- 
placement Sq dq ^ J dqg{q,h)^. In dimensional regularization 
(minimal subtraction scheme) the integral is just the unitarity 
piece, Jo = — 
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for a sharp momentum cutoff A and < p < A, one 
obtains for the phase shifts 



p cot 5 (p) 



MnCo 



2A 

TT 



P 



P , A 

TT A — p 



(26) 



At zero energy, r(0, 0) 
of Co is given by 



2ao/M7r and, thus, the running 



MACo(A) = -- 



ao 



(27) 



2A 



In this case, the phase shift is given by 

pcotd{p) = H — log- 

Q' TT A — p 

= -- 

ao A 



(28) 



= in 



which corresponds to an ERE with = V2 = 
the hmit A ^ oo. Note that finite cutoff corrections scale 
as 1/A. This indicates a relatively slow convergence to- 
wards the infinite cutoff limit and hence that finite cutoff 
effects are quantitatively important and might even be- 
come a parameter of the theory. Actually, one might 
determine A by fixing the effective range from the first of 
Eq. ([28]), ro(A) = 4/(7rA) = 2.77fm, yielding the accurate 
numerical value A = 90.7 MeV. In this case, this par- 
ticular three-momentum regularization method becomes 
itself a model, since we have no control on the remainder. 
In any case, it is straightforward to check that for any 
finite cutoff there is no off-shellness: T{k' , k) = T{p,p). 

The running given by Eq. ([27]) must be used for any 
cutoff A if we want to renormalize in the end. However, 
thinking of the more general case where finite range cor- 
rections are relevant such a running is only reliable for 
very small cutoffs A <C 7r/2ao- If we consider also a C2(A) 
coefficient in the potential, the corresponding LSE can be 
solved with the ansatz 

T{k', k) = n{p) + r2(p)(fc2 + k'^) + T4{p)k^k'^ , (29) 

which yields a set of three linear equations for Tq{p), 
T2{p) and T4^{p). After some algebraic manipulation, the 
final result for the phase shift can then be written in the 
form 



pcot5{p) 



10(C2MA3-K3)V(Af7r) 



9(C|MA5 
2A 



P 1 A 
-log — 

TT A 



5Co) - 15C2(C2MA3 + 6)p2 
P 



■P 



(30) 



Matching at low energies to the ERE, Eq. (fTTI) . we get 
the running of Co and C2 



1 

ao 
1 



10(C2MA3 - 3)2 2A 
9M7r(-C|MA5 -I- 5Co) ~ V 
5OC2 (3 4- C2MA3)^ (6 + CaAfA^ 
277r (-5Co + C^A^Mf 



2 

ttA 



(31) 



The first equation allows to eliminate uniquely Co in 
favour of ao and C2 , but as we see there are two branches 
for the solutions. However, we choose the branch for 
which C2 decouples in the infrared domain, i.e. fulfills 
C2 for A ^ 0. In fact, at small cutoffs, one gets for 
this branch 



AfCo(A)A 
MC2(A)A-'' 



2aoA / 2aoA 



2 ^2aoAy 



(32) 



The factor 2/3 appearing in the small cutoff expansion for 
Co differs already from the coefficient in the case C2 = 0. 
Eliminating Co and C2 in favour of ao and ro, the phase 
shift becomes 



p cot S(j)) = — 



2A 



(tt - 2Aao)2 



Trao 2A(7r - 2Aao) -I- aop2(ro7rA - 4) 



2A p A -hp 
— + - log 

TT TT A — p 

+Irop' + Od). 
ao 2 A 



(33) 



Note that the finite cutoff corrections are, after fixing tq, 
again 0{A~^). So, fixing more low energy constants in 
the contact theory does not necessarily imply a stronger 
short distance insensitivity, as one might have naturally 
expected in other words, the inclusion of a higher 
dimensional operator such as C2 does not improve the 
ultraviolet limit, at least in the polynomial representa- 
tion given by Eq. ()13|) . In Sec. II VI we will show, however, 
that with just one counterterm Co the inclusion of pion 
exchange long distance contributions generates a much 
faster convergence towards the renormalized limit as an- 
ticipated in Refs. [H, 0, \^ (see also Ref. [Ill for a 
quantitative estimate). In Sec.|V]we will also show that 
when a C2 counterterm is added this scaling behaviour is 
not only broken but also the phase shift fails to converge 
in the limit A ^ 00. 

Thus, we see that one can establish a one-to-one map- 
ping between the counterterms Co, C2 and the threshold 
parameters ao and tq. Nevertheless, this is done at the 
expense of operator mixing, i.e., both Co and C2 are in- 
tertwined to determine both the scattering length and 
the effective range. In other words the cutoff dependence 
of Co is different depending on the presence of C2 . As we 
have seen this is not a problem since for small cutoffs we 
expect the running of Co to be fully independent of C2 



We have in mind dispersion relations where any subtraction at 
zero energy and derivatives thereof of the dispersive part improve 
the high energy behaviour and become more insensitive in the 
ultraviolet. As we see this is not the case in the contact pion-less 
theory. 
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and hence on tq. However, unlike the one counterterm 
case, where C2 — 0, the solutions of Eq. ([5^ may become 
complex when 

a^roTrA^ - IGa^A^ + 12ao7rA - Stt^ < . (34) 

For the physical threshold parameters this happens 
already for A > Ac = 382MeV (the other two roots are 
complex). Above this critical value the potential violates 
self-adjointness. For rg ^ one has A^ 16/(7rro) — > 
00. Thus, the cutoff can only be fully removed with a 
self-adjoint short distance potential if vq — 0. This is 
consistent with the violation of the Wigner causality con- 
dition reported in [sl, [13, [H, [HI] . Note that the violation 
of self-adjointness is very peculiar since once Co and C2 
have been eliminated the phase-shift ((33)) remains real 

One feature in the theory with two counterterms Co 
and C2 is that the off-shell T— matrix becomes on-shell 
only in the infinite cutoff limit, 

T{k',k)^T{p,p) + 0{A-') . (35) 

This is unlike the theory with one counterterm Cq where 
there is no off-shellness at any cutoff. Thus, finite cutoff 
effects are also a measure of the off-shellness in this par- 
ticular problem. This will have important consequences 
in Sec. fVl when attempting to extend the theory with two 
counterterms in the presence of the long distance pion 
exchange potentials since the off-shellness of the short 
distance contribution of the potential becomes an issue 
in the limit A — > 00. 

The situation changes qualitatively when the fourth 
order corrections depending on two counterterms C4 and 
C4 are considered. Obviously, we cannot fix both C4 
and C4 simultaneously by fixing the slope parameter V2 
of the effective range expansion, Eq. (fTT|) . Clearly, one 
expects some parameter redundancy between C4 and 
or else an inconsistency would arise since a sixth order 
parameter in the effective range expansion should be 
fixed. The situation worsens if higher orders in the mo- 
mentum expansion are considered due to a rapid pro- 
liferation of counterterms while there is only one more 
threshold parameter for each additional order in the ex- 
pansion. This required parameter redundancy is actually 
a necessary condition for consistency which is manifestly 
fulfilled within dimensional regularization but not in the 
three-momentum cutoff method Moreover, it was re- 
alized some time ago [13, HSj] that the finite cutoff regu- 
larization and dimensional regularizarion in the minimal 



Nonetheless, off-shell unitarity deduced from sandwiching the re- 
lation T-T^ = -2iTiT^5{E-Ho)T between off-shell momentum 
states, is violated, since the Schwartz's reflection principle fails 
[T{E + iO+)] ^ T(E - iO+). This would also have far reaching 
consequences for the three body problem, since three body uni- 
tarity rests on two-body off-shell unitarity and self-adjointness 
of three body forces. 

This operator redundancy has also been discussed on a La- 
grangean level [70l | based on equations of motion and in the ab- 
sence of long distance interactions ( see also iTll V 



subtraction scheme yielded different renormalized ampli- 
tudes for a truncated potential. This non-uniqueness in 
the result due to a different regularization happens when 
a non-vanishing C2 counterterm is considered. In any 
case, the dimensional regularization scheme has never 
been extended to include the long range part of the TPE 
potential which usually appear in the present NN con- 
text. Thus, for the momentum space cutoffs which have 
been implemented in practice the short distance repre- 
sentation is somewhat inconsistent at least for a finite 
value of the cutoff A. 

Alternatively, one may choose an energy dependent 
representation of the short distance physics as 

Vs^Co + 2p'C2 + P^(2C4 + C^) + . . . (36) 

In this case the correspondence between counterterms 
and threshold parameters ao,''o, ^2, etc. is exactly one- 
to-one, and the parameter redundancy is manifest, since 
the on shell T-matrix depends only on the on-shell po- 
tential. Actually, under dimensional regularization the 
representations of the potential Eq. p3)) and Eq. pG]) 
yield the same scattering amplitude. Although this on- 
shell equivalence is certainly desirable it is also unnatu- 
ral, if the long distance potential is energy independent. 
We will nevertheless analyze such a situation in the next 
subsection in coordinate space. 

The previous discussion highlights the kind of undesir- 
able but inherent off-shell ambiguities which arise when 
finite range corrections are included in the short distance 
potential . In our view these are unphysical ambiguities 
which have nothing to do with the unambiguous off-shell 
dependence of the long distance potential. Of course, 
one way to get rid of the ambiguities is to take the limit 
A — > 00 which corresponds to the case where a truly zero 
range theory is approached. However, even for a finite 
cutoff there is a case where one is free from the ambi- 
guities, namely when the short distance potential is both 
energy and momentum independent for s-wave scattering 

Vs{k\k) ^ Co{K) . (37) 

The key point is that we allow only this counterterm to be 
cutoff dependent and real, as required by self-adjointness. 
Of course, the discussion above for the contact theory 
suggests the benefits of using just one Co counterterm 
but does not exactly provide a proof that one must take 
further counterterms such as C2 to zero. The extension 
of this analysis to the case of singular chiral potentials in 
Sec. |V] will yield the definite conclusion that renormaliz- 
ability is indeed equivalent to take C2 — 0. 



^ This fact becomes more puzzling if the potential V = (72 (fc^ -l- 
k'^ — 2p^) is considered. It vanishes on the mass shell k = k' = p 
but nonetheless generates non trivial on shell scattering for the 
three-momentum cutoff. 
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B. Coordinate space 



yielding 



The previous renormalization scheme is the momen- 
tum space version corresponding to the coordinate space 
renormahzation adopted in a previous work by two of us 
(MPV and ERA) [la,!!!!!!. ActuaUy, in the pure con- 
tact theory, we can relate the renormalization constant 
with the momentum space wave function explicitly. At 
large values of the short distance cutoff Vc, the zero en- 
ergy wave function reads, 



uo{rc) = 1 ■ 

ao 



Thus, the following relation holds 



l-rc 



u'oirc) 



ao-rc uoivc) 
Comparing with Eq. (|27)) . we get 



MACo(A) 



u'oira) 



uoirc) 



(38) 



(39) 



(40) 



where the momentum cutoff A and the short distance 
cutoff Tc are related by the equation 



Arr, = 



(41) 



which is nothing but an uncertainty principle relation be- 
tween cutoffs Note that for the standard regular solu- 
tion UQ{r) ^ r one has a vanishing counterterm Co = 0. 
In contrast, Co 7^ for the irregular solution. In the case 
of the singular attractive potentials the solution is reg- 
ular but highly oscillatory and the Co takes all possible 
values for Tc ^ 0. A more detailed discussion on these 
issues can be seen in Refs. '43', '4?, 7^. Of course, strictly 
speaking both Eq. (|40)) and Eq. (|4T|) are based on a zero 
energy state, and in the finite energy case we will assume 
these relations having the limit Tc ^ or A — > 00 in 
mind. 

Let us now deal with finite energy scattering states. 
Since there is no potential, ULir) = 0, for r > Tc we have 
the free wave solution 



Up{r) 



sin(pr + 5{p)) 
sin(5(p) 



(42) 



In the theory with one counterterm, we fix the scattering 
length ao by using the zero energy wave function and 
matching at r = Tc so we get 



uo{rc) Up{rc) 



(43) 



This relation will be shown to hold also in the presence of a local 
potential, see Appendix 1X1 



= pcot{prc + 6{p)) , 



(44) 



and thus 



p cot S (p) 



1 

ao 



1 -pjap - re)tan(prc) 
p{ao - rc) + tan(prc) 

0(re) . 



(45) 



This is in qualitative agreement with the momentum 
space result when the cutoff is being removed, Eq. ([^5]) 
and, as we can see, the approach to the renormalized 
value is similar if the identification rc — tt/(2A) is made. 
Note further that since the boundary condition is energy 
independent the problem is self-adjoint and hence orthog- 
onality between different energy states is guaranteed. 

The theory with two counterterms where both ao and 
To are fixed to their experimental values opens up a new 
possibility, already envisaged in Ref. related to the 
non-uniqueness of the result both for a finite cutoff as 
well as for the renormalized phase-shift. As pointed out 
above, this non- uniqueness was noted first in momen- 
tum space Ref. [stI [s^ when using a finite three dimen- 
sional cutoff or dimensional regularizarion (minimal sub- 
traction). Remarkably, within the boundary condition 
regularization we will be able to identify both cases as 
different short distance representations. 

Actually, when fixing ao and tq wc are led to 



Up{rc) 



dpirc) = 



u'o{rc)+P u'^jrc) 
uo{rc) +p'^U2{rc) 



Oip'). (46) 



where uq and U2 are defined in Sec. Ill CI Note that now 
self-adjointness is violated from the beginning due to the 
energy dependence of the boudary condition. Within the 
second order approximation in the energy the neglected 
terms are 0{p'^), so any representation compatible to this 
order might in principle be considered as equally suitable. 
The close similarity to a Fade approximant suggests to 
compare the following three possibilities for illustration 
purposes 

^ Woirc) + P^u'^jrc) 
^ uoiXc) +p^U2{rc) ' 



jiii 



uoirc) [uo{rc) uo(»'c)^ 

n'oircf ^ 

uo{rc)u'g(rc) +p^ [u2{rc)u'o{rc) - uo{rc)u'2{rc)] ' 

(47) 



and study what happens as the cutoff is removed, 

0. Note that all three cases possess by construction the 

same scattering length ao and effective range tq and no 
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potential for r > r^. Straightforward calculation yields 

pcoi5{p) = -— + Irop^ + O{rl) (I) 
ao 2 



p cot S (p) 
p cot 6 (p) 



— + irop2 + 0(r,) (II) 
"0 2 



1 



1 



ao 1 - ^aor-op^ 



+ 0(re) (III) 



(48) 



As we see, all three representations provide the same 
threshold parameters, but do not yield identical renor- 
malized amplitude for finite energy. Actually, cases I and 
II coincide with the three-dimensional cutoff regulariza- 
tion method (see Sec. lIII Al) . whereas case III corresponds 
to dimensional regularization (MS). Moreover, the finite 
cutoff corrections to the renormalized result are, gener- 
ally, 0{rc) while the rational representation yields correc- 
tions 0{r^). These observations survive at higher orders 
when V2, W3, etc. threshold parameters are further taken 
into account. This indicates that not all short distance 
representations are equally "soft" in the UV-cutoff. The 
generalization of these results to the case of singular TPE 
chiral potentials was studied in Ref. [52| and will be also 
re- analyzed in Sec. IVl while discussing the consistency of 
the standard Weinberg's power counting. 

In any case, when finite range corrections are consid- 
ered within the boundary condition regularization, there 
are two possible renormalized solutions depending on the 
particular parameterization of short distance physics. A 
nice feature of this regularization is that they can be iden- 
tified with similar results found already in the momentum 
space analysis of Ref. [stI [sst when confronting three- 
momentum cut-off and dimensional regularization. We 
note also here that no ambiguity arises when the bound- 
ary condition is assumed to be energy independent, in 
which case self-adjointness is guaranteed. 



IV. RENORMALIZATION OF PION 
EXCHANGES WITH ONE COUNTERTERM 

The study of the contact theory in Sec. |llT] provides 
suggestive arguments why it is highly desirable to carry 
out a regularization with a single counterterm in the ^Sq 
channel by adjusting it to the physical scattering length 
for any cutoff value. In this section we want to extend 
that study when the long distance chiral potential orga- 
nized according to Weinberg power counting enters the 
game and the cutoff is removed. By taking the cutoff to 
infinity, we are actually assuming that all degrees of free- 
dom not included in the present calculation become in- 
finitely heavy. This way we expect to learn about missing 
physics in a model and regularization independent fash- 
ion. The traditional strategy of adjusting an increasing 
number of counterterms may obscure the analysis. In 
other words, by using this minimal number of countert- 
erms, we try not to mock up what might be still missing 
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FIG. 1: LO, NLO, N2LO and N3LO running of the coun- 
terterm (in GeV~^ as a function of the cutoff A in the ^So 
channel for small cutoffs A < 600MeV. The renormalization 
condition is determined by fixing the scattering length to its 
experimental value ao ~ — 23.74fm. We use the parameters 
of Ref. for the pion exchange potential Vl ■ 



in the long range description. 

In this context there is of course the question of con- 
vergence or cutoff insensitivity of the phase shift, when 
A — > 00, provided we keep at any rate the scattering 
length ao to its physical value by suitably adjusting the 
unique counterterm Co (A). In coordinate space this is 
a rather trivial matter if the long distance potential is 
local [H, m, H^, as it happens in the LO, NLO, and 
N2L0 Weinberg counting. The analysis in momentum 
space involves detailed large momenta behaviour of the 
Lippmann-Schwinger equation and one must resort to 
trial and error. Indeed, the N3L0 case analyzed below 
includes nonlocalities and it turns out to provide conver- 
gent results. 

For numerical calculations, we take the values for the 
Ci and di parameters appearing in the pion exchange po- 
tential used in Ref. [2J], which do a good job for pe- 
ripheral waves, where re-scattering effects are suppressed 
and where one is, thus, rather insensitive to cutoff ef- 
fects. We will only consider TPE contributions to the 
N3L0 potential. 



A. Renormalized N3LO-TPE 

To determine the running of the counterterm we start 
from low cutoffs A <^ since the long range part of 
the potential is suppressed and Eq. (|7f)) may be used. 
Actually, the analytical result is well reproduced by the 
numerical method used to solve the Lippmann-Schwinger 
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FIG. 2: LO, NLO, N2L0, and N3LO convergence of the phase shifts as a function of the momentum cutoff A for fixed LAB 
energies, Tlab = 15,60, 105,300 MeV. The renormahzation counterterm Co (A) is always determined by fixing the scattering 
length to its experimental value ao = — 23.74fm. 



equation in the pure contact theory with Cq. Once this 
identification has been done, the value of the cutoff is 
increased steadily so that the scattering length is always 
fixed to the experimental value, ctQ — — 23.74fm. This 
adiabatic switching on of the long range physics guar- 
antees that we are always sitting on the correct branch 
which smoothly goes into the contact theory at low cut- 
offs 

The running of the counterterm Co (A) is depicted for 
the LO, NLO, N2L0, and N3L0 potentials in Fig. [1] in 
the low cutoff region. As expected, the deviations from 
the simple result of the pure contact theory, Eq. (P7)) . 
start at A ~ m,r for LO because of the Itt exchange po- 
tential. For higher cutoffs the NLO, N2L0, and N3L0 
counterterm displays a cycle structure very similar to 
what has been observed in coordinate space [47| . 

The convergence of the ^5*0 phase shift for fixed val- 
ues of the Lab energy, Tlab = 15, 60, 105, 300 MeV is 
displayed in Fig. [2l Of course, one observes a faster con- 
vergence for small energies. For the maximal value of 
Tlab = 300 MeV, cutoff values A - 2 GeV are needed to 
change the phase shift by less than 1°. 

In agreement with the analytical estimates of Ref . , 



^"^ In general, there may appear many solutions for Co(A) fitting 
ctQ. They are physically unacceptable unless they behave as 
MCo(A) — > 2ao/i' for A — > since they do not evolve into the 
theory where the long range components are decoupled. 



the convergence of the regulated phase-shifts towards 
their renormalized values follows a computable power like 
pattern, 5{k) - 6{k,A) = 0(A-"/2-i). The more singu- 
lar the potential at large momenta the faster the con- 
vergence. Thus, the expected increased insensitivity at 
short distances is indeed confirmed. 

Finally, the renormalized ^So phase shift is presented 
in Fig. [3] for LO, NLO, N2L0, and N3L0. As a check of 
the present calculation in momentum space, let us men- 
tion that we re pro duce the coordinate space renormalized 
phase shifts of 'sF] at LO, NLO, and N2L0. For instance, 
at the maximal CM momentum oi p = 400 MeV, the 
maximal discrepancy between the coordinate space and 
momentum space phase shifts is less than half a degree 
when rc = 0.1 fm and A = 4 GeV, respectively. 

The clear converging pattern can be observed all over 
the elastic scattering region, actually at Tlab = 300MeV, 
one has (5^° = 72.72° , (J^^O ^ 6.44°, ^^^lo ^ _24.20°, 
and (5^3LO _ —28.20°. However, there is still a discrep- 
ancy with the Nijmegen PWA result which at this energy 
is = -4.68 ± 0.55° for np scattering. 

To have an idea on the uncertainty of the calcu- 
lated phase shift, we vary the scattering length ao = 
—23.74(2), the value of, gT^NN = 13.1(1) and gA = 
1.26(1). Actually the error in the latter is correlated 
through the Goldberger-Treiman relation, so we will 
quote both as an error in gA only. In addition, for 
the chiral constants we take the central values used in 
[23 | which provided a good description for the peripheral 



12 



> 
O 

II 

< 
to 



80 



60 



40 



20 



-40 









"""<, \ 








OPE 




NLO 




N2L0 




N3L0 




NIJ 





50 100 150 200 250 300 350 400 

p (MeV) 



FIG. 3: LO, NLO, N2LO, and N3LO renormalized phase 
shifts in the '^So channel as a function of the CM momen- 
tum compared to the Nijmegen Partial Wave Analysis p9l |. 
Only one counterterm is used and fixed by the physical value 
of the scattering length. 



waves and for the uncertainties we assume as an educated 
guess those from the TriV study [t^I- The only exception 
is C4, for which the error from t^N [73| is much smaller 
than the systematic discrepancy with the NN determina- 
tion [23|. So it is more realistic to take the systematic 
discrepancy as the error. Thus we take c\ — —0.81(15), 
C2 = 2.80(23), C3 = -3.20 ± 1.35, C4 = 4.40 ± 1.0, 
dx+di = 3.06(21),d3 = -3.27(73), di4-rfi5 = -5.65(41). 
The results for the particular variations are presented 
in table H] and as we clearly see, the uncertainty in C3 
dominates the total error. At Tlab = 300 MeV, one has 
^LO ^ 75.90±0.2°, (JNLO = 6.5±0.7°, J^aLO ^ _24±6", 
^N3LO _ _2g -|- 90, As we see the statistical uncertainty 
stemming from the input parameters increases with the 
order. If we take the difference ^^^lo _ ^N2L0 _ 40 
an estimate of the systematic error, adding in quadrature 
statistical and systematic errors we have, AJ^Q^p^ ~ 10°, 
still a smaller quantity than the discrepancy to the Ni- 
jmegen phase shift. This may suggest that after renor- 
malization there are some physical effects missing even 
at the N3L0 level beyond pure TPE. 



B. Inclusion of Delta 

In the previous section we have seen that all TPE ef- 
fects included to N3L0 display a convergent pattern after 
renormalization, but there is still some missing physics. 
Note that while at LO and NLO the only parameters are 
gA, TOtt, Mat and /tt, at N2L0 and N3L0 there appear 
new low energy constants (the Ci and di, respectively) 



which can be related to ttN scattering and encode short 
range physics not considered explicitly. The A resonance 
is an outstanding feature of TriV scattering and explains 
a great deal of the low energy constants. Thus, it is in- 
teresting to analyze the role of explicit A excitations as 
intermediate states in the NN potential. The importance 
of explicit A degrees of freedom has been emphasized on 
power counting grounds in several previous works with 
finite cutoffs where the A^A split ting is regarded as a 
small parameter ~ 0, [TgI . |27| . |74| |. The crucial role 
played in the renormalization problem has been stressed 
in Ref. [1^. In this section, we analyze the NN poten- 
tial with the NLO terms of Q together with the lA and 
2 A in the box diagrams as computed in Ref. Q. One 
advantage of such an approach is that, as compared to 
the standard A-less theory, there only appears the A^A 
splitting as a parameter. 

The renormalization of the ^6*0 channel proceeds along 
the lines discussed in the A-less theory. The results for 
LO, NLO, NLO-l-lA, and NL0-h2A renormahzed phase 
shifts in the ^Sq channel are plotted and compared to 
the Nijmegen Partial Wave Analysis [i^ as a function of 
the CM momentum in Fig. U) The most striking result is 
the very strong resemblance between the N2L0 A-less vs 
NLO+IA and the N3L0 A-less vs NL0-h2A. These re- 
sults sustain the treatment of the A^A splitting as a small 
parameter ^ 7x1^^ corresponding to the A-counting NLOa 
= NLO + 1A+ 2A. The fact that all the contributions 
fall off at large distances as ^ Q-2m„r gyggggts that de- 
spite the ability to mimic higher order corrections in the 
Weinberg counting for the long distance potential in the 
^5*0 channel there is still missing shorter range physics 
beyond TPE. 



C. Three-pion exchange contributions 

At N3L0, three-pion exchange occurs for the first 
time. These contributions have been calculated by 
Kaiser [l^, |Tl[ and found to be small, which is why 
present N3L0 NN potentials omit these contributions 
when renormalization is not implemented. However, it 
should be noted that, for small distances, the 37r dia- 
grams are proportional to and, thus, will ultimately 
dominate at short distances. A rough estimate of the re- 
sults published in Refs. [13, [HI suggests that the sum of 
all 3tt graphs is attractive. Thus, we infer from there that 
the scattering length will still be a free input parameter. 
On the other hand the 37r-exchange contribution falls off 
as ~ g-3ni„r g^-j- iQj^g distances so it becomes active at 
rather short distances, and so the effect is expected to 
be small because of the short distance suppression (mod- 
ulo oscillations) of the wave function u{r) ~ r''/'* typical 
of potentials with a short distance power like singular- 
ity [H, 13, \^ . This agrees with the rule that the more 
singular the potential the more convergent is the calcula- 
tion, as we have extensively discussed above. Therefore, 
in a complete and renormalized N3L0 calculation of the 
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TABLE L Induced errors in the ^So phase shift (in degrees) at Tlab ~ 300 when the input parameters are varied. The sign 
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^Sq phase shifts, we expect the Stt effects not to be large, 
although the predictions may be closer to the empiri- 
cal values than in Fig. [3l An accurate investigation of 
the impact of Stt exchange at N3L0 on NN phase shifts 
represents an interesting and challenging project for the 
future. 



D. Irrelevance of a fixed C2 counterterm 



As we have mentioned, the very definition of the po- 
tential is ambiguous as it requires fixing the polynomial 
terms in the momentum. In the calculations presented 
above, we have taken the renormalization scheme for the 
potential where the fixed and cut-off independent choice 
C2 = for the potential is made. We have also analyzed 
the situation when a different renormalization scheme is 
taken, namely a non-vanishing arbitrary C2 coefficient 
which does not run with A. However, we allow Co (A) to 
run in a way that the scattering length ao is reproduced. 
This generates a different renormalization trajectory for 
Co (A) as compared to the case C2 = 0. We find by actual 
calculations that this fixed C2 coefficient is irrelevant, 
i.e. the renormalized phase shift does not depend on this 
fixed value in the limit A — s- 00. Roughly speaking, the 
reason is that while the polynomial combination C2q^ is 
large at large momenta, the pion exchange part behaves 
as log (g^) with an additional logarithm and thus dom- 
inates for fixed C2 . This irrelevance of C2 was highlighted 
in the coordinate space analysis of Ref . [3^ where the reg- 
ularization based on a radial cut-off rc would provide a 
compact support not sensing the details of the distribu- 
tional contributions in Eq. ([H]), regardless on how many 
derivatives of delta's are included. A quite different situ- 
ation arises when C2 depends on the cut-off in a way that 
the effective range is fixed as we discuss in the next Sec- 
tion |Vl There, it will be shown that if C2(A) is relevant 
then the phase shift is not convergent. 
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FIG. 4: LO, NLO, NLO+IA, and NLO-H2A renormalized 
phase shifts in the ^ So channel as a function of the CM 
momentum compared to the Nijmegen Partial Wave Anal- 
ysis [49I ]. Only one counterterm is used and fixed by the 
physical value of the scattering length. 



V. WEINBERG'S COUNTING AND 
RENORMALIZATION WITH TWO 
COUNTERTERMS. 

A. Momentum space 

In previous sections, we have seen that the LO, NLO, 
N2L0, and N3L0 chiral potentials can be renormalized 
when one counterterm Co (A) is determined for any value 
of the cutoff A by fixing the scattering length to its ex- 
perimental value, Qfo = — 23.74fm and the limit A —^ 00 
is subsequently taken. This agrees with the observation 
of Refs. [H, 111] that attractive singular potentials 
can be renormalized with a single counterterm. However, 
Weinberg counting requires further counterterms in the 
short distance potential, see Eq. (fT3|) . For instance, both 
at NLO and N2L0 a C2 counterterm should be included. 

In this section, we discuss whether the ^Sq scattering 
amplitude is renormalizable, i.e. whether the scatter- 
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ing amplitude has a well defined limit, when both the 
Co (A) and C2(A) counterterms are included and fixed 
by fitting the scattering length aQ = —23.74 fm as well 
as the effective range tq — 2.77 fm for any value of the 
cutoff A and the limit A oo is pursued. Actually, 
when the NLO (N2L0) long distance potential, Eq. 
is used, this way of proceeding corresponds to renormal- 
izing the NLO (N2L0) approximation of the ^5*0 channel 
in the standard Weinberg counting. The method we fol- 
low in practice is a straightforward extension of the case 
with only just one counterterm Cq. In the case when 
C2 = 0, we get a finite (renormalized) effective range 
r^^'-' = 2.29fm, which is close but still differs significantly 
from the experimental value. At N2L0 with C2 = one 
gets r^2LO _ 2.86 fm. This last value is so close to the 
experimental one that one would not expect big changes 
when a C2 is added to exactly fit the experimental vq. 
Thus, it makes sense to investigate what would happen 
if one uses a non-vanishing C2 to account for the miss- 
ing, and in principle tiny, contribution to the effective 
range. The surprissing result, to be discussed below in 
detail, is that trying to fit the discrepancy in the effec- 
tive range with a C2 counterterm is incompatible with 
renormalizability. 

In the contact theory we found that for A > 380MeV 
the counterterms Cq and C2 diverge before becoming 
complex numbers. Thus, the corresponding Hamiltonian 
becomes non-self- adjoint, in harmony with the Wi gne r 
causality bound violations unveiled in Ref. [5^, ITsj . 
When the NLO contribution to the potential is included, 
the critical value of the cutoff is slightly shifted towards 
the higher values Ac — 480 — 550MeV; above those values 
the short distance contribution to the potential becomes 
non-self-adjoint and causality bounds are violated. Nev- 
ertheless, and similarly to the pion-less theory, the phase 
shifts remain real beyond this critical cutoff. In Fig.[5]the 
value of the phase shift for both NLO and N2L0 approx- 
imations is depicted for several cutoff values. As we see, 
at NLO one observes large variations when the cutoff is 
changed from A = 2000MeV to A = 2780MeV. In Fig. H 
we plot the value at a fixed LAB energies as a function 
of the cutoff for both NLO and N2L0. As we see the 
phase shift does not seem to converge to any particular 
value when A is increased; in fact large variations can be 
clearly seen at moderate cutoff values, A ~ 200MeV for 
Elab = 105MeV and A - SOOMeV for Ei^ab = SOOMeV. 
On the other hand, let us note that a plateau region not 
always appears, and when it does the residual contri- 
bution from pion-exchange effects is less important and 
the counterterms may dominate the calculation. In other 
words, if the cutoff was too small, one would be driven 
back to an effective range expansion with no visible con- 
tribution from chiral potentials whatsoever. Thus, any fi- 
nite cutoff calculation within a higher cutoff regime turns 
out to be strongly cutoff dependent, or else the cutoff 
must be fine tuned to intermediate energy data, hence 
becoming an esssential, and not an auxiliary, parameter 
of the theory. Note that all the problems are triggered 
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FIG. 5: NLO (upper pannel) and N2LO (lower pannel) con- 
vergence of the phase shifts as a function of the CM momen- 
tum for some fixed values of the cutoff A. The renormalization 
counterterms Co (A) and C2(A) are always determined by fix- 
ing the scattering length ao and the effective range ro to their 
experimental values qq ~ —23.74 fm and ro = 2.77 fm. For 
A > Ac ~ SOOMeV both Co (A) and C2(A) become complex 
while the phase shifts remain real. 



by insisting on fitting the effective range parameter to 
the experimental value by introducing the C2 countert- 
erm required by Weinberg power counting on the short 
distance interaction. In contrast, if C2 = 0, as advocated 
in Refs. [H, 0, \^ and Sec. |TV] above, not only is the 
phase shift convergent in the limit A — > cxi (in practice 
A > IGeV) but also most of the effective range is satu- 
rated by the chiral potential 

Although we have checked that these results are nu- 
merically robust within the LSE by increasing the num- 



Actually, in the limit ag ^ 00 and under the assumption of Van 
ver Waals dominance one gets at N2LO the analytical result (see 
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ber of grid points, there is always a reasonable doubt, 
since these are demanding calculations. Thus, it would 
be nice to understand the lack of convergence from a dif- 
ferent perspective, as we do next. 



B. Coordinate space 

Let us compare the previous findings in momentum 
space with related investigations in coordinate space [52[ . 
A thorough study in coordinate space has revealed that 
the very existence of the A — > oo limit may actually de- 
pend on the specific representation of the short distance 
physics. This analytical result has been verified by nu- 
merical calculations and will become extremely helpful 
in analyzing the momentum space calculations presented 
above. Thus, for the sake of completeness it is worth 
reviewing the emerging pattern from Ref. [s^ j. 

As we have discussed in Sec. Ill CI within the bound- 
ary condition regularization the unknown short distance 
interaction is represented by the logarithmic derivative 
of the wave function at the short distance cutoff radius 
Tc, u'f.{rc)/uk{rc). In an energy expansion of the wave 
function at short distances, one has = Uq + fc^U2 -f . . . 
and its logarithmic derivative for which a continuity con- 
dition is required. This introduces an energy dependence 
which eludes the Wigner causality bound discussed in 
Refs. [H, [7^, since self-adjointness is broken from the 
start. As we discussed previously in the contact theory 
described in Sec. IIII Bl at second order in the energy the 
neglected terms are 0{p'^) and any of the three represen- 
tations displayed by Eq. ([iT)) might be equally accept- 
able. Actually, we found that renormalized amplitudes 
fall into two classes and that not all of them approach 
the renormalized limit in the same way, see Eq. (j48p . 
Here, we will extend that study by inquiring what hap- 
pens when the same three representations displayed by 
Eq. (|47)) are used and the long distance potential 
is taken to be the NLO and N2L0 for r > as the cut- 
off is removed, Tc — > 0. Note that all three cases possess 
by construction the same scattering length ao and effec- 
tive range rp and the same long distance potential for 
r > Tc- Thus, any difference is clearly attributable to the 
different short distance representation. The results are 
displayed in Fig. [5] for fixed short distance cutoff values 
Tc as a function of the CM momentum p. For finite values 
of Tc we see a difference which can naturally be explained 
by the different off-shell behaviour of the short distance 
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Ref. [sgl for details), 
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(-4 -I- I5g\ + 24c:iM - 8C4M) 



which yields ro ^ 2.33fm, a surprissingly good approximation. 
Further studies exploiting the chiral Van der Waals correlations 
can be found in Ref. [3511. 
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FIG. 6; NLO and N2L0 convergence of the phase shifts as 
a function of cutoff A for fixed values of the LAB energy, 
Tlab = 105 MeV (upper pannel) and Tlab = 300 MeV 
(lower pannel). The renormalization counterterms Co (A) and 
C2(A) are always determined by fixing the scattering length 
Qfo and the effective range ro to their experimental values 
ao = -23.74fm and ro = 2.77fm. For A > Ac ~ 500MeV 
both Co (A) and C2 (A) become complex while the phase shifts 
remain real. 



physics. As we see, the difference persists as the cutoff is 
being removed and in fact is magnified. As was already 
pointed out in a previous work [5^ , only the case / rep- 
resenting a rational function turns out to yield a unique 
and well defined finite value for the phase shift 

This fact is clearly seen from inspection of Fig. [7] where 
we plot the phase shifts for a fixed value of the CM mo- 
mentum p = 300 MeV as a function of the equivalent 



Note that in the contact theory, this rational representation pro- 
vides the softest regulator, i.e. cutoff effects scale quadratically 



and not linearly as in all others, see Eq. II47I I. 
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FIG. 7: NLO and N2LO running of the phase shifts as a function of the equivalent momentum cutoff A = 7r/(2rc) (in MeV) 
with Tc the short distance cutoff, for fixed values of the LAB energy Tlab = 300 MeV and for several parameterizations of the 
short distance energy dependent logarithmic derivative (see Eq. (|47p . In all three cases, the chiral potential for r > Vc is the 
same and the scattering length qq and the effective range ro are both fixed to their experimental values ao = — 23.74 fm and 
ro = 2.77 fm. 





FIG. 8; NLO (left panel) and N2LO (right panel) convergence of the phase shifts as a function of the CM momentum for fixed 
values of the short distance cutoff Vc for several parameterizations of the short distance energy dependent logarithmic derivative 
(see Eq. (|47p . In all three cases, the chiral potential for r > rc is the same and the scattering length ao and the effective range 
ro are both fixed to their experimental values ao — —23.74 fm and ro — 2.77 fm. 
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sharp momentum cutoff A = n/ (2rc) derived for the pion- 
less theory in Sec. lIIII The striking similarity between the 
momentum space calculation presented in Fig. [6] and the 
coordinate space calculation displayed in Fig. [7] for finite 
cutojfs is noteworthy although not completely surprising 
in the light of the analysis of Appendix El There, the 
finite cutoff equivalence, A = 7r/(2rc), deduced in the 
contact theory is shown to hold also in the presence of a 
local potential. 

Of course, all these features depend on the singular 
character of the interaction at short distances and do 
not depend on the specific form of the potential, so we 
expect them to hold also when A degrees of freedom are 
explicitly taken into account. 

Although we cannot prove it analytically in momentum 
space, the remarkable similarity of the coordinate space 
analysis with the present momentum space calculations 
strongly suggests that the standard polynomial represen- 
tation of the short distance potential is incompatible with 
renormalization. Of course, this does not preclude the 
possible existence of a suitable potential representation 
of the short distance interaction, most likely with com- 
plex counterterms, such that both low energy parameters 
can be fixed and the cutoff can at the same time be re- 
moved. We leave such an interesting investigation for the 
future. 



VI. CONCLUSIONS AND OUTLOOK 

In the present paper, we have analyzed the renormal- 
ization of the singlet np phase shift in the ^So chan- 
nel incorporating one- and two-pion exchange effects. 
For the long distance potential, the standard Weinberg 
scheme based on dimensional power counting is adopted 
to N3L0. However, the short distance physics is param- 
eterized in terms of one unique energy and momentum 
independent counterterm whose cutoff dependence is de- 
termined by adjusting the scattering length to its phys- 
ical value. The present analysis is carried out in mo- 
mentum space in a somewhat complementary manner as 
previously done in coordinate space [4^ . Actually, we 
have analyzed and reproduced those results directly in 
the more popular momentum space after the cutoff has 
been effectively removed. In order to stress the equiva- 
lence of both approaches, we have displayed many results 
in parallel. This is not just a stylistic matter of presen- 
tation; besides numerical simplicity, much understanding 
of the renormalization problem has been achieved by a 
direct analysis in coordinate space at least for local po- 
tentials. The present work provides a further example in 
this respect. 

The momentum space framework allows a direct exten- 
sion to include N3L0 contributions which include non- 
local pieces in the long distance chiral potential. The 
main outcome of such a calculation is that the N3L0 
chiral potential induces rather small corrections as com- 
pared to the N2L0 results in all of the elastic scattering 



region. As anticipated previously [sBIl , this happens to be 
so even when the N3L0 potential is more singular than 
the N2L0 one at short distances. Moreover, the analyti- 
cal scaling behaviour for large cutoffs A of the scattering 
amplitude predicted in a previous coordinate space anal- 
ysis [5^ is confirmed qualitatively by the present mo- 
mentum space calculations; increasing the order in the 
expansion of the long distance potential suppresses more 
the finite cutoff dependence. In summary, the scheme 
is convergent, but differs from the expected phase shifts 
obtained from partial wave analyses. An error analysis of 
the results shows that although there are some uncertain- 
ties induced by the input parameters, the corresponding 
error bands are not large enough to be compatible with 
the Partial Wave Analysis of the Nijmegen group. 

The N3L0 approximation to the long distance chiral 
potential is computed within a heavy baryon expansion 
and assuming only explicit nucleon degrees of freedom 
in the NN potential. This implies in particular that 
the A^A splitting is considered a non-small parameter. 
However, this number is about twice the pion mass, so 
it is not clear whether such an assumption is fully jus- 
tified. Therefore, we have analyzed the NN scattering 
problem when explicit A intermediate state excitations 
are included in the potential. The net result is that, after 
renormalization, the lA and 2A box diagrams contribu- 
tions mimic extremely accurately the A-less N2L0 and 
N3L0 potentials, respectively, and thus fail to describe 
the higher energy region of the ^Sq phase shift. Given 
the fact that these contributions also fall off at large dis- 
tances as Q-'^rn^r ^ ^j^jg j-gg^j^ rciuforces the conclusion 
that there is some shorter range missing physics beyond 
that provided by TPE. 

Further, it is interesting to compare the present renor- 
malized results with those obtained for a two pion ex- 
change potential computed within a rclativistic baryon 
framework which sums up all heavy nucleon compo- 
nents [ll] thus contains a full determination of TPE con- 
tributions assuming all other degrees of freedom (includ- 
ing explicit A's) are infinitely heavy. There, all np partial 
waves with j < 5 are analyzed and for the particular case 
of the ^5*0 channel the results are rather similar to those 
found here. This, again, suppports the present conclusion 
that, presumably, the components of the NN potential 
with a shorter range than TPE might finally provide the 
missing repulsion needed to reduce the 3% overshooting 
of the ^5*0 effective range as well as the too low value of 
the phase shift in the region with CM momenta p > m^r . 

We have also investigated the convergence of the stan- 
dard Weinberg counting for both the short distance as 
well as the long distance potential. At NLO and N2L0 
this corresponds to including a polynomial momentum 
dependence in the short distance interaction by means 
of two counterterms which may then be fixed by adjust- 
ing the scattering length and the effective range to their 
experimental values. We find that, with these two renor- 
malization conditions, the counterterms turn out to be 
complex for not too large cutoffs (around 500MeV), sig- 
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naling the breakdown of self-adjointness of the potential, 
and more generally suggesting a violation of Wigner's 
causality condition. Nonetheless, despite phase shifts re- 
maining real, a unique renormalized limit does not exist. 
This lack of convergence agrees with similar findings in 
coordinate space anticipated in Ref. [E^l- The physical 
understanding of the situation is as follows. On the one 
hand the long range physics is fixed, thus, the off-shell 
propagation within the long range region is unambigu- 
ous. On the other hand, the short range potential must 
be adjusted to provide the threshold parameters such as 
the scattering length, the effective range etc. If the short 
range potential is just determined from fixing the scatter- 
ing length only, this is at zero energy and becomes zero 
range when the cutoff is removed, so there is no off-shell 
ambiguity. In contrast, fixing further the effective range 
requires non-zero energy and there are in fact infinitely 
many ways of parameterizing this, with different off-shell 
behaviour even when the cutoff is removed. In other 
words, fixing a finite range and removing the cutoff de- 
pends on details of the method. These conflicts between 
off-shcllness and finiteness are not new in field theory. 
Green's functions which are renormalized on-shell do not 
necessarily provide finite off-shell amplitudes. The stan- 
dard Weinberg parameterization as a polynomial in mo- 
menta is one choice which may not turn out to be con- 
sistent with renormalizability. In this regard, it is inter- 
esting to mention that coordinate space detailed studies 
based on renormalization group properties (s^ suggest 
the existence of a suitable short distance representation 
yielding convergent results in theories with more than 
one counterterm. The generalization of those results to 
momentum space would require an in-depth study of the 
renormalization group in the presence of eventually non- 
local but singular potentials. We note that the energy 
dependence of the coordinate space solution violates self- 
adjointness explicitly and we expect that most likely an 
energy independent non-polynomial momentum space so- 
lution would spoil self-adjointness as well. 

It is important to note that with a finite cutoff and four 
counterterms the ^5*0 phase shift has been successfully de- 
scribed within the standard Weinberg counting 24, 25]. 
In this regard, it is natural to question the usefulness of 
taking the infinite cutoff limit and to carry out a renor- 
malization process. From a physical point of view the 
limit A ^ oo corresponds to consider other degrees of 
freedom than those not considered explicitly to be in- 
finitely heavy. Thus, the aim of the renormalization 
program is far more stringent than previous finite cut- 
off calculations. Within such a framework any failure 
can be unambiguously attributed to missing physical in- 
formation on the long distance potential, and the renor- 
malization process highlights this in a rather vivid man- 
ner as one clearly sees when going from OPE to TPE 
potentials. From a mathematical perspective the renor- 
malizability requirement imposes tight constraints on the 
admissible forms of the short distance physics for some 
given long distance interactions. These are powerful con- 



ditions which have traditionally been the great strength 
of the renormalization ideas to avoid unnecessary prolif- 
eration of interactions. Since the power counting of long 
distance potentials is not uniquely determined yet (one 
may e.g. include or not explicit A's), it would be very 
helpful to see what kind of short distance constraints on 
those potentials might be imposed as well on the basis of 
renormalizability or other principles. 

The problems with the Weinberg counting for the short 
distance polynomial form in the momenta of the interac- 
tion with renormalization found in this paper for NLO 
and N2L0 is to be added to the other ones noted in pre- 
vious works [13, H^. On the mathematical side, it is 
noticeable that the coordinate space analysis of Ref. [s^ 
conjectured this result by exploiting the compelling re- 
quirement of completeness and self-adjointness for the 
renormalized quantum mechanical problem for a local 
chiral potential which at first sight may seem completely 
germane concepts to the EFT machinery. On the phe- 
nomenological side, it should be noted that this incon- 
sistency result does not explain why the renormalized 
phase shift with just one counterterm comes out reason- 
ably close to the accepted ones, but certainly makes this 
unique and finite prediction more credible and inevitable 
from a theoretical perspective if self-adjointness is main- 
tained, and provides further confidence on the virtues of 
a renormalization principle within the chiral approach to 
the NN interaction. 
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APPENDIX A: MOMENTUM VS COORDINATE 
SPACE FOR FINITE CUTOFFS AND LOCAL 
POTENTIALS 

In this appendix we discuss further the relation be- 
tween momentum and coordinate space for finite cut- 
offs and local potentials. We will show that solving the 
Lippmann-Schwinger equation for a long range poten- 
tial with a sharp cutoff A is equivalent to the solution 
of the Schrodinger equation for the same potential in a 
discretized grid with A?' = tt/A. This equivalence re- 
minds the previous identification between sharp momen- 
tum cutoff A and the short distance cutoff rc found in the 
contact theory in Sec. Illli A = 7r/(2rc). We can choose 
the initial grid point at Vc = 7r/(2A), thus recovering the 
previously mentioned equivalence. Actually, by invoking 
Nyquist theorem on the cut-off LS equation, we deduce 
a discretized version of the variable phase equation fn\ 
which allows to discuss both the renormalization as well 
as the decimation problem of the NN-force based on chi- 
ral potentials. 



1. The cutoff Lippmann-Schwinger equation 

Let us consider as a starting point the Lippmann- 
Schwinger equation for s wave scattering, written in the 
form 

T{k',k) = Vik',k) + M I dqV{k',q) f , Tiq,k), 

Jo p - q 

(Al) 

where T{k' , k) are the matrix elements of the T— matrix 
between initial and final CM momentum states k and 
k' respectively and the corresponding potential matrix 
element is given by 

2 1"°° 

V{k',k) = - / drjo{k'r)jo{kr)V{r)r^ , (A2) 



where d{x) is the Heaviside step function. In coordinate 
space the cutoff potential becomes 



for a local potential. Since we are integrating the in- 
termediate momentum up to infinity, we are implicitly 
assuming that V{k',k) is a regular potential. The sin- 
gular potential case will be discussed later. Now, if we 
cutoff the potential in momentum space we get the 
regularized potential 

VA{k', k) = 9{A - \k'\)6{A - \k\)V{k', k) , (A3) 



Cutting off the high momentum states is not exactly the same 
integrating out the high energy states, which produces a low mo- 
mentum effective energy dependent "optical" potential. We are 
focusing on the long range potential here. The short missing dis- 
tance piece could be included by a Taylor expansion in momenta 
or energies (see Sec. I A i\ . 



^A(r',r) 



k'^dk / k''^dk'jo{kr)jo{kr')V{k,k') 



(A4) 



which is obviously nonlocal and becomes local only when 
A ^ cxo (second line). The cutoff LSE, where \q\ < A, can 
be solved by standard means. In the spirit of an EFT, 
based on the idea that low energy dynamics does not 
depend on the details at short distances, it may perhaps 
be appropriate to proceed a bit differently. Actually, if 
high momentum states are cutoff from the theory the 
same idea should apply to small resolution scales, i.e. 
regardless whether they are short or long. That means 
that using too much information on the potential V{r) 
even if it is exactly known point-wise may be illusory at 
wave-lengths longer than a given resolution, Ar, so 
that one can sample V{r) with some resolution Ar ~ 
1/A. This obviously reduces the number of coordinate 
mesh points in the integration. 



2. Nyquist theorem 

In the present context, Nyquist theorem [tI] is remark- 
ably useful because it provides an optimal way of sam- 
pling signals which have a bandwidth in Fourier space, 
i.e., functions for which 



e"'='f{x)dx = for |fc| < A , (A5) 



then, for the original function we have 



/(x) = / e-^^-F{k)^ 

-A 



(A6) 



If we define the sampling function fs{x) of the function 
f{x) at the equidistant points Xn — n/S.x, with the opti- 
mal Aa; — -K I K separation, 

oo 

fs{x) = ^ j{x^)l\x5[x - Xn) , (A7) 

n— — oo 

and compute its Fourier transform we get 

oo 

Fsik) = Ax J2 /(a;n)e^'""- (A8) 



n— — oo 



For example, if a potential V{r) is supplemented by highly os- 
cillatory ripples at short resolution scales, Ar the phase shifts 
should be insensitive to them at A 2> Ar. 
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Note that for the optimal samphng Fs{k) = F{k) 
Inverting the Fourier transform we get 



fsix) 



-A 



Ztt 



E /(.„) ""[f^"-f] A.. (A9) 



Thus, if we sample the function according to Eq. (IA7p . 
the following identity holds at the sampling points 



fsiXn) = fiXn) ■ 



(AlO) 



Hence, there is no loss of information on the sampling 
points if the sampling is done equidistantly with the op- 
timal Nyquist frequency. Ax = tt/A. In particular, it 
does not really make sense to sample the function with 
smaller Ax. In the next section we apply this sampling 
principle to the potential. 



3. The optimal grid for sampling the potential in 
coordinate space 

Nyquist theorem also applies for the special case of 
the Lippmann-Schwinger equation in s-wave scattering, 
although the derivation is different in some details to the 
one presented in the previous section. For clarity, we 
present here these details. 

We will consider first the general case of a non-local 
potential, for which the Schrodinger equation in s-wave 
reads 

/■OO 

-u"{r) + M / dr V{r,r)u{r') ^ k^u{r) . {All) 
Jq 

We can sample this non-local potential as 

OO 

Vs{r,r') - {Arf ^ l^(r„„ r„) <5(r - r„0 <5(r - r„) . 

m,n— 

(A12) 

For the momentum space representation of the sampled 
potential we get 

„ OO 

Vs{k,k') = {Arf- V V^.njo{kr^)jo{k'rn), 

IT ^ ^ 



m,n— 



(A13) 



As the original F(k) is bandwidth limited to the [—A, A] interval, 
it can be expressed as a Fourier sum 



nk)= E 



a„ e 



ikxji 



where Xn = titt/A. Thus it is trivial to see that Fs{k) = F{k) 
when the sampling is done with Ax = tt/A. 



where Vm,n = rmTn V{rm,Tn)- As a consequence of the 
cut-off A, the potential V{k, k') can be expressed as a 
sum of spherical bessel functions 

oc 

v{k,k')^ E 

^n.m Jo 

(fcr„)jo(fcV„), (A14) 

m,n— 

where r„ = titt/A. Then by taking Ar = tt/A, the sam- 
pled potential recovers the original one, i.e. Vs{k,k') = 
V{k, k'), for fc, k' < A. By Fourier-transforming Vs{k, k') 
back to coordinate space, it can be checked that it repro- 
duces the original sampling points, i.e. 

Vs{rn,rm) = (^Ar) V(r„, r„) = y(r„, r„) (A15) 

for the Nyquist sampling frequency. 

In the case of a local potential, the one which interests 
us most, we sample the following way 

OO 

r -rn). (A16) 



After double Fourier-transforming, we get 



A 



Vs(?'„,r„) = F(r„)— (5, 



(A17) 



which makes the potential local for the grid points. It is 
in fact a finite cut-off version of V{r, r') = V{r) 6{r — r') 
once we notice that — 6nm 5(r — r'). 



4. The Lippmann-Schwinger equation with a 
cut-off 



For a Lippmann Schwinger equation with a finite cut- 
off A all the matrix elements become a linear combination 
of separable terms. Thus the LS equation becomes a 
linear matrix equation which can be solved by standard 
techniques by writing 



T{k,k') = ^Ti.ijo{kn)jQ(krj). 
V{k,k') = E^*J-J"(^^^)jo(^'^j) 



(A18) 
(A19) 



and defining propagator matrix elements 



Gy= / dq- jjo('?r,)jo(grj)- (A20) 



Therefore we get 

Tij — Vij + E ^uGimT, 



ni] 



(A21) 



Im 



This equation can be reduced to a finite-dimensional N x 
N linear algebra problem by cutting the sums to i, j — N. 
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The effect of this siniphfication can be seen by samphng 
Vij in coordinate space (see Eqs. (|A7p and (|A16p ) 



N 



Vir) ^ Ar^\^(r„)<5(r-r„) = — ,(A22) 



n=0 



from which we can check that cutting the sum is equiv- 
alent to introducing the (harmless) infrared cut-off r^r. 
Note that for local potentials, the sampled potential ma- 
trix elements at the grid points are diagonal 



Thus Eq. (|A2ip becomes 



(A23) 



(A24) 



which looks like a multiple scattering equation, with on- 
shell propagation between delta-shell scatterers. After 
matrix inversion, the on-shell solution is then given by 

T{p) = J2 T^MUpri)k{prj) ■ (A25) 



5. The discretized Schrodinger equation 

In the previous section we have seen that for the 
Lippmann-Schwinger equation with a cut-off A, we 
can either use the original momentum space potential 
V{k^ k') and solve by standard means, or expand this 
potential, i.e. use the sampled potential Vs{k,k')^ and 
solve as a linear algebra problem (as they are both the 
same potential for momenta below the cut-off). 

Alternatively we can directly solve the Schrodinger 
equation for the sampling potential Vs{r). This proce- 
dure will give an excellent approximation to the solu- 
tion of the LS equation (although not the exact solution, 
as explained at the end of this section) but at a much 
smaller computational cost. For this purpose we make 
the replacement 

N 

V{r) ^ Vs{r) = ^r^r - n)V{n) , (A26) 

i=0 

which is a superposition of equally spaced delta-shell po- 
tentials. It should be noted that the point ro = does 
not contribute, as it lies at the integration boundary of 
the Schrodinger equation. Then, we have 



Ar 



TT 

A 



TN =i NAr. 



(A27) 



Thus, for a potential of size a where we do not want to 
describe energies higher than A we should do with an 
infrared cut-off larger than the potential's length scale. 



tn ^ a, or equivalently A'' ^> Aa/ir. We can solve the 
Schrodinger equation piecewise. 



u{r) = Aism{kr + 5i^i/2) 
n-i <r<ri, 



(A28) 



where Ai is the amplitude and (5j„i/2 can be understood 
as the accumulated phase-shift due to adding a new delta 
shell at Vi chosen to be located at the midpoint r^^i 
(for reasons to become clear soon). Note that with this 
choice the lowest possible location of the phase-shift, S1/2 
corresponds to the point 



TT 

2A 



(A29) 



which we may identify with a short distance (ultraviolet) 
cutoff. This identification between the momentum-space 
cutoff agrees with the one obtained for the pure short 
range theory by solving the LSE without any discretiza- 
tion (see Sec. IIIip . 

It should be noted that the previous method for solv- 
ing the Schrodinger equation does not really generate the 
exact solution of the LS equation with a cut-off, but only 
a close approximation. Although Vs{r) is the best sam- 
pling function for V\{k,k'), it is not the same quantum 
mechanical potential. The solution to the LS equation 
with V\{k,k') can be exactly reproduced by solving the 
non-local Schrodinger equation, Eq. (jAlip . with the po- 
tential Vs{r, r') which comes from inverse Fourier trans- 
forming Vs{k^k'). The Schrodinger equation should be 
solved with trivial initial conditions, u(0) = 0, as all the 
physically relevant information is included in V5(r, r'). 

On the contrary, if we solve Vs{r), a sum of delta shells, 
we must include a non trivial initial condition even 
in the absence of any short range physics. If we solve 
the discretized Schrodinger equation with a trivial initial 
boundary condition, 5i/2 — 0, then a certain error will be 
included in the final solution. In the worst case we can 
expect an error of order Ar. By means of the variable 
phase equation [t^ we can perform a better assessment 
of the error yielding, for example, an ©((Ar)^) error 
for a Yukawa potential or ©((Ar)"^) for a square well. 

From the previous discussion, it is apparent how to 
include a counterterm in the computation. It enters 
throught the initial sampling point of Vs^ which maps 



Specially since we are effectively ignoring the ro = sampling 
point. 

For small enough Ar, the variable phase equation will yield 

5i ~ -k / M V(R) dR . 
2 Jo 

As a curiosity, we can see that 5\/2 scale as an inverse power 
of Ar for a singular potential, thus signaling the need of a coun- 
terterm. 
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onto a constant term in Vs{k, k') Thus it is ignored 
when solving the Schrodinger equation for Vg, and it en- 
ters through the initial condition 5i/2. Similar remarks 
can be made for singular potentials, in which a countert- 
erm must be included in order to obtain a stable result 
for Ar 0. 



6. The discrete variable phase equations 

The previous discussion can be elaborated further to 
reach interesting results. Matching the wave functions at 
the points where the delta shells are located, r = r.i, we 
simply get 

kcot{kn + S,+i/2)-kcot{kn + di_i/2) = ArUin). 

(A30) 

where U{r) — 2iJ,V{r) ~ MV{r) is the reduced poten- 
tial. This is a recursion relation for the phase-shift at 
the interval midpoint r.^^i — {i + i)7r/A. Unlike the ma- 
trix equation which has traditional storage limitations 
for large number of grid points N, this equation does not 
posses this shortcoming allowing for rather large N val- 
ues. Similar equations were deduced many years ago [t^ 
as a practical tool to attack the inverse scattering prob- 
lem and to determine the NN potential on the grid points. 
Defining the discretized effective range function, 



Mi = k cot Si 



we get 

M,- 1 1 fc cot kri — k 
Af,-+ 1 + k cot kri 



M.ifccot kri — k 

' 2 

M, _ 1 + fccot kri 



(A31) 

= ArU{r,) 
(A32) 



which can be rewritten as a continuous fraction. Note 
that for the cutoff theory both Eqs. (jA30P and (IA32j) are 
almost exact. The only approximation comes from the 

For example, the Cq counterterm when projected to the s-wave 
takes the form 



Co 



Co 



(when discretizing). 



Thus, if we write the sampled potential Vs(k,k') 



6 a, 



0,i+i 



3r, 



6 (ao,^+i - n) 



finiteness of the cutoff A. An important property which 
will be used later on is the reflection property, namely 
the symmetry under the replacement 



Ar 



-Ar 



(A33) 



which means that on the grid running the relation up- 
wards or downwards are inverse operations of each other. 
Obviously this property may fail in practice only due to 
accumulation of computer round-off errors over large evo- 
lution distances. This is also the reason why we choose 
the midpoint 7'i+i/2 for the accumulated phase shift: in 
any other case we would lose the reflection property. 
The regular solution at the origin reads 



S^{k)=Q S^^Uk) = S{k) 



(A34) 



If we take the limit A ^ oo we can define S{k,ri) — (5,;(fc), 
to get 



dS{k,R) 
dR 



1 



U{R) sin^ {kR + S{k,R)) 



0{Ar^). 
(A35) 

which is the variable phase equation [77[ up to finite 
grid corrections and can be interpreted as the change 
in the accumulated phase when a truncated potential of 
the parametric form U{r)6{R — r) is steadily switched 
on as a function of the variable R. Eq. (|A30p is thus a 
discretized variable phase equation, corresponding to a 
discretized potential sampled at the optimal Nyquist fre- 
quency. This equation and its generalization to coupled 
channels has extensively been used to treat the renormal- 
ization problem in NN scattering in Refs. (isl. [t^. 

The discrete equations for the low energy parameters 
can be easily computed from the low energy expansion 



k cot 5i{k) = — + ^roik'^ 



- V2Ak* 



(A36) 



which yields 



yZjj Vij jo(kri)jo(krj), and take into account Eq. | |A23| |. we see 

that Co maps only onto Vbo = ^~^> giving a constant contribu- 
tion in momentum space, as expected. 



ArU{r,) (A37) 
6ao,i-i - + a2 J^_6r, + 3ro ^_ 



(A38) 



and similar equations for higher low energy parameters, i'2,i'3,i'4, • • ■ • Note that the potential only enters ex- 
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plicitly in the equation involving the discretized scatter- 
ing length. One of the nice features of this equation is 
the way how it handles the case of singular points, when 
oq^R) or other parameters diverge, since standard inte- 
gration methods for the corresponding differential equa- 
tion are based on the smoothness of the solution and 
hence fail. 

The way to proceed in practice is quite simple. For 
a given number of grid points N we take the scattering 
length, Qq = «0j the effective range, r-Q^jv -1.1/2 — 

and run Eq. (|A37|) and Eq. (jA38|) to determine ao,i/2 = 
Q!o(7r/2A) and ro,i/2 ~ ro{n/2A). Then, we use, e.g. 
Ml/2 = ~1/q^o,i/2 +^0.1/2^^/2 and run Eq. (|A32[) up- 
wards to obtain M ^ k cot S = M, 



N+l/2- 



Due to 



the reversibility of the algorithm one exactly has M 
— 1/ao + rofc^/2 in the limit fc — > for any finite grid 
size Ar (modulo computer arithmetic round-off errors). 
In this way we can fix the initial boundary conditions to 
exactly reproduce any given scattering length, effective 
range, etc. 

This method was used successfully [13] to extract low 
energy threshold parameters in all partial waves with j < 
5 from high quality potentials [5l| . On that case, we ran 
Eq. I|A37P and Eq. (|A38|) with trivial initial boundary 
conditions ao^i/2 = 0, ^0.1/2 0, etc, and obtain the 
potential's threshold parameters as ao = o.o,n+i/2i ''o = 
''o,Ar-i-i/2, etc. 

As we have said the discretized running parameters lo- 
cated at the lowest possible radius Vc = 7r/2A correspond 
to the short distance interactions. Nowhere in the equa- 
tions does the potential at the origin appears explicitly, 
since according to Eq. (jA32[) . one starts with [/(tt/A). 
From this point of view the treatment of regular and sin- 
gular potentials at the origin is on equal footing. 



7. The decimation process in momentum space 



imation problem of the LSE where the low energy states 
are cutoff. Assuming that we have the LSE with a given 
cutoff A related to potential, we ask how does the physi- 
cal phase-shift change when we make the transformation 
A — > A/2, in the potential fixed. Applying the process ex- 
plained above based on Nyquist theorem we see that this 
corresponds to double the grid resolution Ar — > 2Ar. 
Obviously, by repeating the process we may effectively 
have Ar ^ a (being a the range of the potential) and 
hence the potential never contributes and i5Ar+i/2 = '^i/2- 
Thus, for Aa/n ^ 1, the short range theory is recovered. 
In the opposite limit Ao/tt ^ 1 we get the full short 
range plus long range physics. 

If we now set a grid with a fixed number of points 
NAr = Ntt/K ~ a, then making A A/2 corresponds 
to a — > 2a so it could be viewed as switching on the 
potential. 

The discretized representation of the LSE is alge- 
braically closed, since it corresponds to a separable inter- 
action, but any term contains all powers in momentum. 
Moreover, they also vanish for large k, unlike any poly- 
nomial approximation to it. The numerical procedure 
explained above of computing the phase-shift from some 
low energy parameters like a,ro, etc. and the long dis- 
tance potential V{r) can be also carried out by assuming 
that the potential at the lowest grid point, Vc = tt/2A, is 
energy dependent yielding 



Vs{k,k') = Vip) JO 



Trk 
2A 



Jo 



vrfc' 
2A 



(A39) 



Expanding in powers of momentum wc can identify the 
terms 



Vsik,k')^V{p) 



12A2 



{e + k'') 



(A40) 



The previous equations provide the accumulated phase 
shifts due to the addition of equidistant delta-shells sam- 
pling the original potential in coordinate space. We want 
to show that they actually provide a solution of the dec- 



with Co, C2, etc. Note that the radius of convergence 
of the expansion in k and k' is the whole complex plane, 
due to the meromorphic character of the spherical Bessel 
functions. 
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